o 

(N 



O 
(N 






a 



o 



> 

rn 

in 

o 
o 



X 



5 pages, 3 figures 

Spatial distribution of local tunneling conductivity due to interference and Coulomb 
interaction effects for deep and shallow impurities on semiconductor surfaces 

V. N. Mantsevicli:] and N. S. MaslovfO 
Moscow State University, Department of Physics, 119991 Moscow, Russia 

(Dated: May 21, 2010) 

Spatial distribution of local tunneling conductivity was investigated for deep and shallow impuri- 
ties on semiconductor surfaces. Non-equilibrium Coulomb interaction and interference effects were 
taken into account and analyzed theoretically with the help of Keldysh formalism. Two models were 
investigated: mean field self-consistent approach for shallow impurity state and Hubbard-I model 
for deep impurity state. We have found that not only above the impurity but also at the distances 
comparable to the lattice period both effects interference between direct and resonant tunneling 
channels and on-site Coulomb repulsion of localized electrons strongly modifies form of tunneling 
conductivity measured by the scanning tunneling microscopy/spectroscopy (STM/STS) . 
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I. INTRODUCTION 



Impurity states at surface and interfaces of semicon- 
ductors strongly modify the local electronic structure 
and consequently determine the behavior of tunneling 
characteristics in STM/STS contacts [l|-|3|. Experimen- 
tal and theoretical investigations of tunneling through 
impurity atom energy level in the case of multichannel 
transport reveal Fano-type line shape in local tunneling 
conductivity when STM metallic tip is positioned above 
the impurity [J, la, |8| . Transformation of Fano-type line 
shape in local tunneling conductivity depending on dis- 
tance value from impurity was also recently investigated 
[3, 01 . Most of the experiments are carried out with 
the help of scanning tunneling microscopy/spectroscopy 
technique [1, 8 10] and theoretical calculations usually 
deals with Green's functions formalism [3, lli|. Com- 
parison between experimental STS results obtained at 
different distance values from the impurity [5| and the- 
oretical investigations of single impurities influence on 
tunneling conductivity [7| provide information weather 
electron transport occurs coherently or incoherently and 
gives an opportunity to initialize impurity type. All these 
effects are caused by local changes of the initial density of 
states due to interactions of non-equilibrium particles in 
the contact area. Taking into account Coulomb interac- 
tion of conduction electrons with non-equilibrium local- 
ized charges can result in nontrivial behavior of tunneling 
characteristics calculated in the case of STM metallic tip 
positioned above the impurity atom [Ill4l3l | . However in- 
fluence of Coulomb interaction effects on local tunneling 
conductivity measured apart from the impurity is a prob- 
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lem of great interest due to the possibility of impurities 
types initialization. 

So in this work we present the modification of formula, 
which describe spatial distribution of local tunneling con- 
ductivity in vicinity of impurity in the case of interference 
between resonant and direct tunneling channels [7| due 
to Coulomb interaction. We performed calculations for 
two extreme cases when resonant tunneling takes place 
through deep impurity state and through shallow im- 
purity state. We applied mean field self-consistent ap- 
proach for shallow impurity state and in the case of 
deep impurity non-equilibrium Coulomb interaction ef- 
fects were studied with the use of Hubbard-I model. Both 
approaches were analysed with the use of Keldysh for- 
malism [IJ]. We have found that taking into account 
Coulomb interaction in addition to interference between 
tunneling channels leads to drastical changing of the tun- 
neling conductivity form depending on the values of tun- 
neling rates and distance from impurity. 



II. THE SUGGESTED MODEL AND MAIN 
RESULTS 

We shall analyze tunneling between semiconductor sur- 
face {ID atomic chain) and metallic STM tip for deep 
and shallow impurities in the presence of Coulomb in- 
teraction. The model of tunneling contact formed by 
semiconductor and metallic tip is depicted in Fig. [TJ IZ? 
atomic chain consists of similar atoms with energy levels 
El and similar tunneling transfer amplitudes 3 between 
the atoms. Distance between the atoms in the chain is 
the same and equal to a. Atomic chain includes impu- 
rity with energy Ed, tunneling transitions from the im- 
purity atom to the semiconductor and metallic tip are 
described by the tunneling transfer amplitudes t and T 
correspondingly. Direct tunneling between the surface 
continuous spectrum states and tip states is described 



by the transfer amplitude t. 
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FIG. 1: Schematic diagram of direct and resonant tunneling 
channels. 



The model system semiconductor-impurity state- 
metallic tip in the presence of Coulomb interaction can 
be described by the Hamiltonian: H: 



H — Hq + Himp + Htun + H^ 



tip 



kc kk' <y 

k p kp 

Himp = 2J ^dC'^a'^da + Uudand^a', Htip = ^J ^P^pa^^iiD 
da p<y 



Indices k and p label the states in the left (semiconduc- 
tor) and right (tip) lead, respectively. The index d in- 
dicates that impurity electron operator is involved. U 
is the on-site Coulomb repulsion of localized electrons, 
nda = c^^Cda, Cda dcstroys impurity electron with spin 
a. Hq is a typical Hamiltonian for atomic chain with 
hoppings without any impurities. Htun describes reso- 
nant tunneling transitions from the impurity state to the 
semiconductor and metallic tip and direct transitions be- 
tween the tunneling contact leads. Himp corresponds to 
the electrons in the localized state formed by the im- 
purity atom, Htip describes conduction electrons in the 
metallic tip. 

With the use of diagram technique for non-equilibrium 
processes [1J| one can get expression for the spatial distri- 
bution of local tunneling conductivity without Coulomb 
interaction [7|: 



dl 



\J lkplkdlpdvlReG^^{uj)cos{2k,j: (cj).t) + 
Ikpilkd + jpd)i^prnGfj^{uj)cos{2kx{u>)x) 



+ 



Ikdlpd 



■ImGU^) + 



Ikd + Ipd 

lkdlpdlkp{i^ - ed)cos{2kx{uj)x) 






+ 



£d)^ + ilkd + IpdW 
(w - ed)2 + 72^ - 7pd(w 



£d) 



{lo - EdY + {"fkd + IpdY 
{uj - Ed)^ + i'jkd + 7pd)^(l - cos{2kj;{uj)x)) 



(w - Ed)^ + ilkd + Ipd)"^ 

Ipdijpd + jkd)cos{2kx{uj)x) , 



(w - £d)2 + {-Jkd + IpdY 



-) 



(2) 



where impurity retarded Green's function is defined by 
the expression: 



G 



dd 



W 



i^ - £d - iilkd + Ipd) 



(3) 



Relaxation rates jkd, Ipd are determined by electron tun- 
neling transitions from localized states to the leads k and 
p continuum states and relaxation rate ^kp corresponds 
to direct tunneling transitions between k and p contin- 
uum states EpT^I^^Gpp = Ipd', Ept^I^Gpp = Ikp^ 

j:kr'/rnGl^k=nd. 

v^ is unperturbed density of states in semiconductor. 
Expression for kx{u}) can be found from the ID atomic 
chain dispersion law: 



i^ikx) = 25 • cos(A:a:a) 



(4) 



Let us now analyze modification of formula [2] due to 
Coulomb interaction effects in the case of deep impuri- 
ties with the help of Hubbard-I model (Hamiltonian of 
suggested model has the form of expression [T] ) . Taking 
into account Coulomb interaction in this model leads to 
formation of two well separeted impurity energy levels Ed 
and Ed + U instead of one initial level Ed- It is reasonable 
to use approximation in which the strongest interaction 
of the considered model - the on-site Coulomb repulsion 
U - is included in G'^^°'{u!). So we can write down ex- 
pression for Gj^'^{uj): 






1 



w 



I](c.) = 



nd- 



LU - Ed- S(w) 
jU{uj-Ed) 



Ed- [l -nd-a)U -i6 



(5) 



With the help of Keldysh diagram technique [Tj] using 
system of Dyson equations one can get expression for 
impurity retarded Green's function. 



GB.cr 
dd — 



1 



UJ -Ed- S(w) - i{jkd + Ipd) 



(6) 
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FIG. 2: Local tunneling conductivity as a function of applied bias voltage calculated for different values of the distance x 
from the deep impurity atom along the atomic chain in the absence (gray line) and in the presence (black line) of Coulomb 
interaction U. For all the figures values of the parameters a = 1, t = 0,22, T = 0,18, t — 0,50, SJ = 1,00, Ed = —1,00, 
U — 3, 50, u?. = 1 are the same. 



To evaluate the dependence of local tunneling conduc- 
tivity on the distance from the deep impurity in the pres- 
ence of Coulomb interaction expressions for ImG^j^(Ld) 
and ReG^^{uj) in formula [H have to be retarded by ex- 
pressions for ImGfj{uj) and ReG^^{uj). It is also nec- 
essary to take into account the presence of electron spin, 
it results in extra summation over a and leads to ad- 
ditional factor equal to 2 in formula [5] Calculating ex- 
pressions for ImG^J (uj) and ReG^^{uj) requires solution 
of self-consistent system of equations. Two of them are 
equations for ImG^^ '^{^)j calculated from expression 
ini and three equations determine impurity atom non- 
equilibrium electron filling numbers: 



nd^o- = / du>nd^cri^)ImG^^'^{uj,nd±cr) 

, . , . Ikdnliuj) + -/pdu^piu}) 
nda[^) = nd-a[i^) = ; 

Jkd + Ipd 



(7) 



energy levels positions is determined by the value of a 
distance. 

It is clearly evident that in the case of deep impurity 
state Coulomb interaction effects lead to negligible shift 
of a resonant peculiarity which corresponds to value of 
applied bias voltage equal to initial impurity energy level 
position (w = Sd)- Suggested model also demonstrates 
absence of resonant peculiarities spreading with increas- 
ing of relaxation rates and distance value from the impu- 
rity. 

Now let us analyze modification of formula [2] due to 
Coulomb interaction effects in the situation of shallow 
impurity state. In this case Coulomb interaction effects 
can be taken into account with the help of mean-field 
approximation (adopted parameters of the model corre- 
spond to the mixed valence regime). It means simply 
that impurity energy level position depends on Coulomb 
interaction of the non-equilibrium electron density. New 
impurity energy level position improved by Coulomb in- 
teraction can be found from the equation: 



where f^fc(w) and n^(uj) are equilibrium filling numbers 
in the tunneling contact leads. 

Figure [5] shows tunneling conductivity as a function 
of applied bias voltage for different values of distance 
from the deep impurity. Tunneling conductivity calcu- 
lated above the impurity has a resonant dip when ap- 
plied bias voltage is equal to impurity energy level posi- 
tion (oj ~ Ed) both for taking (Fig. [2^ black line) and not 
taking into account Coulomb interaction effects (Fig. [5^ 
gray line). Tunneling conductivity calculated above the 
impurity in the presence of Coulomb interaction also re- 
veals peak at the value of applied bias voltage equal to 
initial impurity energy level position shifted on the value 
of Coulomb potential {uj = Sd + U) in comparison with 
the case when Coulomb interaction effects are neglected. 

At the fixed parameters of tunneling contact existance 
of a dip or a peak in tunneling conductivity in both reso- 
nances when applied bias voltage is equal to the impurity 



£d = £d + U{nd) 



(8) 



where Sd is the initial position of impurity energy level 
without Coulomb interaction. Now the main point is 
that the non-equilibrium electron filling numbers nd for 
impurity atom must satisfy self-consistency condition: 



nd 



— / dujnd{uj)ImG^^ 



(cu) 



(9) 



where ImG^^{uj) is determined from equation |3] con- 
sidering Ed substitution instead of Ed- Impurity filling 
numbers Ud^uj) can be found from kinetic equations for 
Keldysh functions G< 0: 



nd{uj) 



-ykdnljuj) + ■^pdn°{uj) 
Ikd + Ipd 



(10) 
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FIG. 3: Local tunneling conductivity as a function of applied bias voltage calculated for different values of the distance x 
from the shallow impurity atom along the atomic chain in the absence (gray line) and in the presence (black line) of Coulomb 
interaction U. For all the figures values of the parameters a = 1, t = 0,22, T = 0, 18, r = 0,40, Sy = 1,00, Sd = —0,15, 
U — 1, 00, u^l = 1 are the same. 



After calculating impurity atom electron filling num- 
bers Ud from the system of equations [3J [51ITU1 we can de- 
termine new position of impurity energy level e^. Next 
step in local tunneling conductivity calculation in the 
presence of Coulomb interaction for shallow impurities 
consists in replacement of energy level value Ed by new 
value ?d in formula [2] 

Figure [3] shows tunneling conductivity as a function of 
applied bias voltage for different values of distance from 
the shallow impurity. In this case tunneling conductivity 
calculated above the impurity has Fano line shape due to 
interference between resonant and non-resonant tunnel- 
ing channels both for taking (black line) and not taking 
(gray line) into account Coulomb interaction effects (Fig. 
[3^) . Tunneling conductivity calculated without Coulomb 
interaction (Fig. [3^ gray line) shows a resonant dip when 
applied bias voltage is equal to impurity energy level po- 
sition (w = Ed)- Considering Coulomb interaction for 
shallow impurities in the mean field approximation leads 
to a resonant dip shift to the higher values of applied 
bias voltage (Fig. [3^ black line). This effect also takes 
place with increasing of the distance value from the im- 
purity moreover shift is accompanied by dip's spreading 
(Fig. [5)d black line). Dip's spreading is a result of both 
impurity initial energy level movement due to Coulomb 
interaction effects and increasing of relaxation rates in 
the studied system. It was found that dip's spreading 
is mostly significant when relaxation rate value exceeds 
value of initial impurity energy level position Sd- Fur- 
ther increasing of distance (Fig. [St black line) shows 
disappearing of the resonant dip in the case of Coulomb 
interaction opposing to the calculations carried out with- 
out taking into account Coulomb interaction effects (Fig. 



[St gray line) . Disappearing of a dip is a result of nearest 
to the dip peaks spreading (one of the peaks corresponds 
to the lower value of applied bias voltage and another 
to the higher value of applied bias voltage in comparison 
with dip position). 

III. CONCLUSION 



In this work we have analyzed the role of interplay 
effects between interference in tunneling processes and 
Coulomb interaction in spatial distribution of tunneling 
conductivity. We have studied two extreme cases when 
resonant tunneling takes place through deep impurity 
state and through shallow impurity state. In the case 
of deep impurity state non-equilibrium Coulomb inter- 
action effects were analyzed with the help of Hubbard-I 
model. For shallow impurity state Coulomb interaction 
of localized electrons was treated self-consistently in the 
mean field approximation. 

For shallow impurity state taking into account 
Coulomb interaction leads to the shift of a resonant dip 
to higher values of applied bias voltage and to the dip's 
spreading with increasing of tunneling relaxation rates 
values in comparison with the case of Coulomb interac- 
tion neglecting. 

For deep impurity state taking into account Coulomb 
interaction results in formation of additional peculiarity 
which corresponds to the value of applied bias voltage 
equal to initial impurity energy level position shifted on 
the value of Coulomb potential. 

This work was supported by RFBR grants and by the 
National Grants for technical regulation and metrology. 
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